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We suggest using a two-color evanescent light field around a subwavelength-diameter fiber to trap 
and guide atoms. The optical fiber carries a red-detuned light and a blue-detuned light, with both 
modes far from resonance. When both input light fields are circularly polarized, a set of trapping 
minima of the total potential in the transverse plane is formed as a ring around the fiber. This 
design allows confinement of atoms to a cylindrical shell around the fiber. When one or both of 
the input light fields are linearly polarized, the total potential has two local minimum points in the 
transverse plane. This design allows confinement of atoms to two straight lines parallel to the fiber 
axis. Due to the thin thickness of the fiber, we can use far-off-resonance fields with substantially 
differing evanescent decay lengths to produce a net potential with a large depth, a large coherence 
time, and a large trap lifetime. For example, a 0.2-/^m-radius silica fiber carrying 30 mW of 1.06- 
^m- wavelength light and 29 mW of 700-nm- wavelength light, both fields are circularly polarized 
at the input, gives for cesium atoms a trap depth of 2.9 mK, a coherence time of 32 ms, and a 
recoil-heating-limited trap lifetime of 541 s. 

PACS numbers: 32.80.Pj,32.80.Lg,03.75.Be,03.65.Ge 



I. INTRODUCTION 

One of the key problems of matter-wave physics is trap- 
ping and guiding neutral atoms . Atom traps and atom 
waveguides can be used as tools for atom optics, atom 
interferometry, and atom lithography 0. Atoms can be 
trapped and manipulated by the gradient forces of light 
waves Q. In particular, the evanescent light waves have 
been used extensively to trap and guide atoms since they 
have high spatial gradients and use rigid dielectric struc- 
tures likeprisms and fibers to define the potential shape 
0, 0, 0, Q ■ An example is a hollow optical fiber with 
light propagating in the glass and tuned far to blue of 
atomic resonance 0- Inside the fiber, the evanescent 
wave decays exponentially away from the wall, produc- 
ing a repulsive potential which guides atoms along the 
center axis. Alternatively, a red-detuned light in the hol- 
low center of the fiber can also be used to guide atoms @ . 
In several experiments , cold atoms have been trapped 
and guided inside a hollow fiber. 

Recently, a method for trapping and guiding neutral 
atoms outside a thin optical fiber has been proposed 
The scheme is based on the use of a subwavelength- 
diameter silica fiber with a red-detuned light launched 
into it. The light wave decays away from the fiber wall 
and produces an attractive potential for neutral atoms. 
To sustain a stable trapping and guiding outside the fiber, 
the atoms have to be kept away from the fiber wall. This 
can be achieved by a centrifugal potential barrier I3| . 

Another way to produce a trap with a potential mini- 
mum outside a surface is to employ the idea of Ovchin- 
nikov et al. 0, who proposed the use of two col- 
ors (i.e., red and blue detunings) and differing evanes- 
cent decay lengths to obtain both attractive and repul- 
sive forces. The two-color method has been considered 



for planar prisms |^, dielectric microspheres free- 
standing channel waveguides ^ , and integrated optical 
waveguides T^. In these systems, the net trapping po- 
tential is small compared to the optical potentials of the 
red- and blue-detuned lights. Such traps are sensitive to 
small field perturbations. As already pointed out by Bar- 
nett et al. |ll| , such sensitivity would be greatly reduced 
if we could increase the difference between the evanes- 
cent decay lengths of the red- and blue-detuned lights. 
In addition, utilizing large detux iing s would also be ad- 
vantageous for coherent guiding pTlll2|. 

In this paper, we demonstrate that a subwavelength- 
diameter optical fiber carrying a red-detuned light and a 
blue-detuned light can be used to trap and guide atoms 
outside the fiber. Due to the thin thickness of the fiber, 
we can use far-off-resonance lights with substantially dif- 
fering evanescent decay lengths to produce a net poten- 
tial with (1) a deep minimum, (2) a large coherence time, 
and (3) a large trap lifetime. We consider two schemes 
of input field polarization. In the scheme where both 
light fields are circularly polarized at the input, a set of 
trapping minima of the potential in the transverse plane 
is formed as a ring around the fiber, and the atoms can 
be confined to a cylindrical shell around the fiber. In 
the scheme where one or both of the input light fields 
are linearly polarized, the potential has two local mini- 
mum points in the transverse plane, and the atoms can 
be confined along two straight lines parallel to the fiber 
axis. 

Before we proceed, we note that, due to recent de- 
velopments in taper fiber technology, thin fibers can be 
produced with diameters down to 50 nm 0, . There- 
fore, a two-color trap using a subwavelength-diameter 
fiber is a quite realistic design. In addition to the po- 
tential practical applications for trapping and guiding 



of atoms, a tapered fiber with an intense evanescent 
field can also be used as an atomic mirror Gen- 
eration of light with a supercontinuum spectrum in thin 
tapered fibers has been demonstrated The evanes- 
cent waves from zero-mode metal-clad subwavelength- 
diameter waveguides have been used for optical obser- 
vations of single- molecule dynamics Rigorous cal- 
culations for the spatial distribution, the waveguide dis- 
persion, and the polarization orientation of the light field 
around a thin fiber have been reported 0, 0, 0| . Thin 
fiber structures can be used as building blocks in the fu- 
ture atom and photonic micro and nano devices. 

The paper is organized as follows. In Sec.|n]we study 
the scheme where both input light fields are circularly 
polarized. In Sec. IIIII we study the scheme where one or 
both of the input light fields are linearly polarized. Our 
conclusions are given in Sec. IIVI 



II. TWO-COLOR TRAP WITH CIRCULARLY 
POLARIZED INPUT LIGHTS 

Consider a thin single-mode optical fiber that has a 
cylindrical silica core of radius a and refractive index ni 
and an infinite vacuum clad of refractive index 77.2 = 1 . 
Such a fiber can be prepared using taper fiber technology. 
The essence of the technology is to heat and pull a single- 
mode optical fiber to a very thin thickness maintaining 
the taper condition to keep adiabatically the single-mode 
condition 19] . Due to tapering, the original core is al- 
most vanishing. Therefore, the refractive indices that 
determine the guiding properties of the tapered fiber are 
the refractive index of the original silica clad and the 
refractive index of the surrounding vacuum. The refrac- 
tive index and the radius of the tapered silica clad will 
be henceforth referred to simply as the fiber refractive 
index ni and the fiber radius a, respectively. 



A. Two-color optical potential 

Consider an atom moving in a potential U outside 
the fiber. If U has a local minimum at a point out- 
side the fiber, the atom can be trapped in the vicinity 
of this point. Our task is to create a potential with 
a trapping minimum sufficiently far from the fiber sur- 
face to make the effects of surface interaction and heat- 
ing negligible. For this purpose, we use two light fields 
in the fundamental modes 1 and 2 with differing fre- 
quencies LOi and [^2, respectively (wavelengths Ai and 
A2, respectively, and free-space wave numbers fci and 
k^, respectively). A schematic of our design is shown 
in Fig. n Assume that the single-mode condition |23| 
Vi = kia^J n\{LOi) — n\ < Vc = 2.405 is satisfied for both 
modes (j = 1 or 2). In addition, assume that the input 
fields are circularly polarized. In this case, the polar- 
ization of the transverse component of each propagating 
field rotates elliptically in time, the orbit rotates circu- 
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FIG. 1: Scliematic of atom trapping and guiding around an 
optical fiber. 

larly in space, and the spatial distribution of the field in- 
tensity is cylindrically symmetric |l8i |. Outside the fiber, 
in the cylindrical coordinates {r, tp, z}, the time-averaged 
intensity of the electric field in mode i is given by 0, |23 

- £l[Kl{q,r) + w,Kl{q,r) + hKl{q,r)]. (1) 

Here the notation Kn stands for the modified Bessel func- 
tions of the second kind, qi characterizes the decay of 
the field outside the fiber, Wi and ft describe the devia- 
tions of the exact fundamental mode HEn from the ap- 
proximate mode LPqi , and £i determines the strength of 
the electric field. The coefficients wt and fi are defined 
as = 2g2//?2(i _ 5^)2 ^nd /, = (1 + - s,f, 

where Si = {\/qja? + l/h'ja^)/[J[{h.ia)/hiaJi{hia) + 
K[{qia) / qiaKiiqitt)]. The parameters qi and hi are 
related to the longitudinal propagation constant f3i as 
q, = {(3f - nifc2)i/2 and = K(L^,)fc2 _ pf^/^ The 

parameter Pi is determined by the eigenvalue equation 
for the fundamental mode at the frequency uji |20| . 

We assume that the atom is in the ground state and 
the fields are off resonance with the atom. The optical 
potential of the atom in the field of mode i is then given 
by ^ 

U, = ~^a,\E,f, (2) 

where ai — a{uji) is the real part of the atomic polar- 
izability at the optical frequency uJi. The factor 1/4 in 
Eq. (01 results from the fact that the dipole of the atom 
is not a permanent dipole but is induced by the field, giv- 
ing 1/2, and from the fact that the intensity is averaged 
over optical oscillations, giving another 1/2. 

In general, the function a{uj) for a ground-state atom is 
given by M «('^) = (eV^e) faji^l " " 
cj^)^ + Ija^'^]- Here e and uie are the electric charge 
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and mass, respectively, of the electron, and Uja^ faj, 
and are the frequency, absorption oscillator strength, 
and emission transition probability, respectively, of the 
spectral line ja. When we use the relation jja = 
faje'^uj'j^ga/2TTmeeoc^gj, we find 



a{uj) — 27reoc' 



,3 9 J 



7ja(l-wV^?a) 



j 



9a (W -a 



(3) 



Here gj is the statistical weight of the excited level 
and ga is the statistical weight of the ground-state man- 
ifold \a). 

We note that, in the model of a two-level atom, the 
real part of the polarizability can be approximated as 
a — —TTCoc^jba/^jJba^ = —(P/hA. Here A = w — uJba is 
the detuning of the optical frequency uj from the atomic 
frequency LUba and d is the projection of the dipole mo- 
ment onto an axis. In deriving the above approximation 
for a, it has been assumed that A is large compared 
to the linewidth jba but is small compared to the opti- 
cal and atomic frequencies. The corresponding approxi- 
mate expression for the optical potential of the atom is 
Uq = hrt^/A where J7 = d\E\/2ti is the Rabi frequency. 
When the field frequency uj is far from resonance with 
the atomic frequencies u>ja, we must take into account 
the multilevel structure of the atom. 

Assume that the timescale of atomic motion is much 
slower than the beating period of the two light fields, 
that is, the inverse of their frequency difference. Then, 
the optical potentials of the two fields add up to give the 
net optical potential 



U = Ui+U2 



(4) 



The sign of the optical potential of each mode is con- 
trolled by the sign of the mode detuning. We choose 
a red detuning (Ai < 0) for the field in mode 1 and a 
blue detuning (A2 > 0) for the field in mode 2. This 
allows both attractive (red-detuned) and repulsive (blue- 
detuned) potentials to be created. When we substitute 
Eq. (P) into Eq. (j^Jl and then the result into Eq. Q), we 
obtain 

U{r) - G2[K^{q2r)+W2Kf{q2r) + hKi{q2r)] 

- Gi[Kl{q^r) + w^Kliqiv) + fiKHq^r)]. (5) 

Here d = ai£l/A = \ai\£l/A and G2 = -a2f|/4 = 
1 0:2 1 £"1/4 are positive coupling parameters. They are pro- 
portional to the powers of the corresponding light fields. 
Since wi < W2, we have qi < q2, that is, the evanes- 
cent decay length Ai — l/qi of the red-detuned field is 
larger than the evanescent decay length A2 = 1/92 of the 
blue-detuned field 8, 20] . When the intensity of the blue- 
detuned field is large enough, the net optical potential U 
is repulsive at short range and attractive at long range 
0. Such a potential possesses a local minimum point, 
where the two forces cancel each other, leading to the 
possibility of atom trapping. Since U is independent of 



ip and z, a minimum point in the one-dimensional space 
{r} corresponds to a cylindrical shell of minimum points 
in the three-dimensional space {r, (p, z}. 

To get insight into the trapping properties of the net 
optical potential U , we perform an analytical treatment. 
To simplify this treatment, we assume that, from one 
side, the coefficients Wi and fi are sufficiently small and, 
from another side, the parameters qiU are sufficiently 
large that the contributions of the terms containing Wi 
and fi in Eq. jSJ are not substantial Then, we ob- 
tain the following approximate expression: 



Uir) = G2K^oi<l2r) - GiK^oilir)- 



(6) 



We use this approximation only for our analysis but not 
for our numerical calculations. 

We introduce the notation r„i and = U{r„i) for 
the local minimum point and local minimum value of J7, 
respectively. From U'{r„i) — 0, we find 



Firm) = A, 



where 



and 



F(r) = 



Ko{qir)Ki{qir) 
Ko{q2r)Ki{q2r) 



92(^2 



(7) 



(8) 



(9) 



Since qi < 92, the function F{r) is a monotonically in- 
creasing function of r. Hence, Eq. Q has a solution 
rm> a\i A> F{a), that is, if 



G2 ^ qiKo{qia)Ki{qia) 
Gi q2Koiq2a)Ki{q2a)' 



(10) 



Condition (|1U|I means that there exists a trapping mini- 
mum outside the fiber if the intensity of the blue-detuned 
light is large enough or if the intensity of the red-detuned 
light is small enough (but not zero). Note that, when the 
power of the blue-detuned light increases or the power of 
the red-detuned light decreases, the depth Ud = —Um = 
\U,n\ of the trapping potential decreases and the local 
minimum point r„i is shifted farther away from the fiber 
surface. Oppositely, when the power of the blue-detuned 
light decreases or the power of the red-detuned light in- 
creases, the potential depth increases and the local min- 
imum point is shifted towards the fiber surface. 

The atom can be prevented from coming too close to 
the fiber surface before being loaded into the trap if the 
net optical potential achieves a non-negative value at the 
fiber surface, that is, if U{a) > 0. This requirement will 
be fulfilled if 



G2 ^ K^jqia) 
Gi - Kl{q2a)' 



(11) 
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When condition is satisfied, condition H10(l is also 
satisfied. To maximize the depth of the trapping mini- 
mum, we optimize the powers of the two field in such a 
way that condition Hill) reduces to an equality, namely. 



G2 ^ Kliq.a) 
Gi Kl{q2ay 



Then, Eq. ® yields 



U{r) = Go 



K^iq^a) Kliq^a) 



(12) 



(13) 



where Gq = GiKl{qia) = G2K^iq2a). 

We analyze the effects of the evanescent decay param- 
eters qi and q2 on the trapping potential under condi- 
tion (|12|l . For this purpose, we fix qi and Gi (or (72 and 
G2) and increase q2 and G2 (or decrease qi and Gi) in 
such a way that equality H12|l is kept. This leads to a 
reduction in U and, consequently, to a decrease in the 
minimum value Urm that is, to an increase in the depth 
Ud = —Urn of the trapping potential. Such variations 
also make the minimum point of the potential closer to 
the fiber surface. Note that, when the fiber radius is 
small, the evanescent decay length A is a fast increasing 
function of the light wavelength A Therefore, using 
a thin fiber, we can obtain a large ratio of q2 to qi by 
choosing a small A2 and/or a large Ai. Thus, a reduction 
in the blue-detuned light wavelength A2 or an increase 
in the red-detuned light wavelength Ai allows us to in- 
crease the potential depth provided the power of the light 
field E2 or i^i, respectively, can be optimized appropri- 
ately to satisfy condition (|12|l . However, A2 must not be 
too small because it is limited by the single-mode con- 
dition. In addition, Ai must not be too large because a 
larger wavelength Ai leads to a larger penetration length 
Ai and, therefore, to a smaller potential magnitude for 
the same power of light. Furthermore, the detunings 
of the light fields from the atomic resonance frequency 
must not be too large because a larger detuning leads to 
a smaller polarizability and, consequently, to a smaller 
potential magnitude for the same power of light. 

In general, the polarizability of an atom is a complex 
characteristic. The imaginary part of the polarizability 
is given by 



kIoj) = 27reoc^ 7^-/ 2 



(14) 



It is responsible for spontaneous scattering. The rate of 
spontaneous scattering caused by a single light field Ei 
is given by 



r 



(sc) 



1 

in' 



- K,- 1 Ei 



(15) 



where = Spontaneous scattering limits the co- 

herence time of the trap. For atoms spending time close 



to the minimum r„ 
scattering rate is 



of the two-color potential f/, the net 



rsc = r(^''(r™)-i-rf^(r^) (i6) 

and the characteristic coherence time is 



1 



'^coh 



(17) 



(rcc) 



Every scattered photon imparts a recoil energy Q 
{hki)'^/2M to the atom, where M is the mass of the atom. 
Therefore, the absorption of mode photons and emission 
of other photons result in a loss of atoms from the trap- 
ping potential. For a trap depth Ud, the quantity [l2 | 



Ud 



(18) 



characterizes the trap lifetime due to recoil heating. 
When the light field frequencies are near to the atomic 
resonances, the scattering rates are large and, therefore, 
the coherence time and the trap lifetime are small. To 
produce a trap with a large coherence time and a large 
trap lifetime, we must use far-off-resonance fields. 



B. van der Waals potential 

An atom near the surface of a medium undergoes a van 
der Waals force. The van der Waals potential of an atom 
near the surface of a cylindrical dielectric rod is given by 



00 ^00 



V{r) = 



n— — 00 



+ {k^ +n^lr^)Kl{kr)] / dea(^OGn(^0, 

(19) 



where 



G„(w) = 



[e{uj) ~ eo]In{ka)I'^{ka) 



eoInika)K'Jka) - t{uj)r^{ka)Kn{kay 



(20) 



Here e(w) is the dynamical dielectric function and /„ is 
the modified Bessel function of the first kind. 

We calculate the van der Waals potential of a ground- 
state cesium atom near a silica fiber surface. The 
dynamical dielectric function of silica is given by 



^ 1 I 0.6961663 A'' , 0.4079426 , 8974794 A^ 

Co ^ ' ^A^ -0.0684043^ ^ A^-0. 1162414^ ~^ A'-^-9. 896161^ ' 

where A is in the units of /xm. To calculate the inte- 
gral p9(l . we use the approximate expression a{uj) = 
27reoc^Ej5i7ja/[5aw|„(w|„ - tJ^)]. This approxima- 
tion is justified because the resonant frequencies of the 
ground-state cesium atom are substantially different from 
the resonant frequencies of silica. We take into ac- 
count four dominant lines of the atom, namely, Aia = 
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FIG. 2: van der Waals potential V of a ground-state cesium 
atom near a thin cylindrical silica fiber. 
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FIG. 3: Ratio between the van der Waals potentials V and 
Vfiax of a ground-state cesium atom near a thin cylindrical 
silica fiber and a fiat-surface bulk silica dielectric, respectively. 



852.113 nm, Aza = 894.347 nm, Asa = 455.528 nm, 
and A4Q = 459.317 nm, see |2^. The emission transi- 
tion probabilities of these lines are 71a = 3.276 x 10^ 
s"i, = 2.87 X 10^ s-\ 73a = 1.88 X 10^ s-i, and 
74a = 8 X 10^ s~^. The statistical weights of the four 
corresponding upper states are gi = 4, 32 = 2, (73 = 4, 
and 34 = 2. The statistical weight of the ground state is 
<7a =2. 

In Fig. 13 we plot the van der Waals potential of a 
ground-state cesium atom near a cylindrical silica fiber 
as a function of the atom-to-surface distance D — r — a. 
The comparison between the solid line (a — 0.2 /im) and 
the dashed line (a = 0.4 //m) shows that a smaller fiber 
radius a leads to a smaller magnitude and a less steep 
slope of the van der Waals potential V. 

When the fiber radius a is small compared to the 
atom-to-surface distance D, the van der Waals poten- 



tial V is, in general, different from the flat-surface 
bulk- medium van der Waals potential Vflat = —Cz/D^. 
Here the coefficient C3 is determined by C3 — 
(?i/167r^eo) /o°° '^^a(*C)[e(*C)-eo]/[e(jC) + eo]- For cesium 
atoms and flat silica surfaces, this coefficient is estimated 
to be C3 9^ 5.6 X 10-*^ J m^ ^ 4.1 x 10^^ mK ^m^. Fig- 
ure 13 illustrates the difference between the fiber-surface 
potential V and the flat-surface potential Vflat- The fig- 
ure shows that V/Vnat < 1, that is, the magnitude of V is 
smaller than the magnitude of T4iat- When the atom-to- 
surface distance D tends to zero, the ratio V/Vfiat tends 
to unity, that is, the two potentials V and T4iat tend to 
become the same. When D increases, the ratio V/Vnat 
reduces, that is, the relative difference between V and 
Vfint increases. 

C. Total potential 

The total potential C/tot of the atom is the sum of the 
net optical potential U and the van der Waals potential 
V, i.e., 

C/tot ^U + V. (21) 

We use Eqs. ©, lO, and (EU to calculate the total 
potential C/tot of a ground-state cesium atom outside a 
thin cylindrical silica fiber with two circularly polarized 
input light fields. The ground-state cesium atom has two 
strong transitions, at 852 nm {D2 line) and 894 nm {Di 
line). To trap the atom, we use red- and blue-detuned 
lights with wavelengths Ai = 1.06 /im and A2 — 700 
nm, respectively. The detunings of the lights from the 
dominant D2 line of the atom are Ai/27r = —69 THz 
and A2/27r = 76 THz. The fiber radius is a = 0.2 /im. 
This radius is small enough to create a large ratio of q2 
to qi and to reduce the effect of the van der Waals force 
.8]. For the above parameters, we find qia = 0.2438 
and (72Q — 0.9686. The corresponding evanescent decay 
lengths are Ai = 0.8 /im and A2 = 0.2 /im. The relative 
difference between the decay lengths is measured by the 
parameter q;a = (Ai — A2)/A2 = 3. The obtained value 
of this parameter is much larger than the characteristic 
value a A = 0.47 estimated for the case of channel guides 
with the TE and TM modes [111- 

In Fig. 21 we plot the contributions to the total trap- 
ping potential C/tot of the atom. The figure shows that 
the evanescent decay length of the red-detuned light is 
substantially larger than that of the blue-detuned light. 
The net optical potential U has a local minimum value 
Urn — —2.9 mK at = 0.37 /im, well outside the fiber. 
The van der Waals potential is only about —5 /iK at 
the distance — a = 0.17 /im from the fiber surface 
and makes just a small contribution in the vicinity of 
the trapping minimum. It is interesting to note that the 
depth Ud — 2.9 mK of the net optical potential is com- 
parable to the red-detuned potential and is larger than 
the blue-detuned potential at the trapping point. These 
characteristics of the two-color fiber scheme are much 
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FIG. 5: Transverse-plane profile of the total potential Utot 
produced by lights that are circularly polarized at the input. 
All the parameters are the same as for Fig. 2] 



FIG. 4: Contributions to the total trapping potential fAot of 
a ground-state cesium atom outside a vacuum-clad subwave- 
length silica fiber with two circularly polarized input light 
fields. The radius of the fiber is a = 0.2 nm. The absolute 
value of the red-detuned component Ui (dotted line) subtracts 
from the blue-detuned component U2 (dashed line) to give the 
net optical potential U (thin solid line). This is modified by 
the van der Waals surface interaction, giving the total poten- 
tial Utot (thick solid line). The light wavelengths are Ai — 1.06 
fim and A2 = 700 nm. The light powers are Pi = 30 mW and 
P2 = 29 mW. 



better than those of the two-color integrated-waveguide 
scheme, where the red- and blue-detuned potentials are 
large but the net potential is small 0. 

We estimate some critical trapping parameters for the 
case of Fig. 0] We find that the rates of scattering 
due to the trapping fields at the potential minimum 



Accord- 



are = 22.39 s-i and r^''' = 8.46 s 

ingly, the net scattering rate and the coherence time are 
Tsc — 30.85 s~^ and Tcoh — 32 ms, respectively. Since 
the recoil energies due to single red- and blue-detuned 
photons are 9^'^'^'' = 0.06 /iK and 62°'^'' = 0.15 /iK, re- 
spectively, the trapping lifetime due to recoil heating is 
estimated to be rtrap — 541 s. At the trapping minimum, 
the radial oscillation frequency is cj^/Stt = 492 KHz, giv- 
ing an estimate of about 23.6 /iK for the energy of the 
atomic mode spacing. 

Due to the circular polarization of the input fields, J7tot 
is cylindrically symmetric. In Fig. |S1 we plot the spatial 
profile of C/tot in the fiber transverse plane. The figure 
shows that Utot is independent of the azimuthal angle tp 
and has a set of deep minima formed as a ring surround- 
ing the fiber. Since the atom can move freely along the 
direction z of the fiber axis, Utot is also independent of z. 
Therefore, the ring of the potential minima in the trans- 
verse plane results in a cylindrical shell of the minima of 
the trapping potential in the three-dimensional space. 

Figure illustrates the effect of the power Pi of the 
red-detuned light on the total potential Utot ■ As seen, an 
increase in Pi leads to an increase in the depth Ud of the 
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FIG. 6: Effect of the power Pi of the red-detuned light on the 
total potential Utot ■ The power of the blue-detuned light is 
fixed at P2 = 29 mW. All the other parameters are the same 
as for Fig. 0] 



trapping potential and to a shift of the local minimum 
point 7'm towards the fiber surface. However, the increase 
in Pi also reduces the height of the repulsive wall in the 
region of r < r^. 

Figure |7| illustrates the effect of the power P2 of the 
blue-detuned light on the total potential Utot- As seen, 
an increase in P2 leads to a decrease in the depth Ud of 
the trapping potential and to a shift of the local minimum 
point r„i farther away from the fiber surface. Meanwhile, 
the height of the repulsive wall in the region of r < r„j is 
increased. The comparison between Figs. El and [7| show 
that the powers Pi and P2 of the red- and blue-detuned 
lights, respectively, have opposite effects on the total po- 
tential of the atom. Therefore, to produce a potential 
with a deep trapping minimum outside the fiber and with 
a high repulsive wall in the region of r < r^, the ratio 
between Pi and P2 must be optimized appropriately. 

Due to the cylindrical symmetry of the total poten- 
tial Utot, the component of the angular momentum 
of the atom is conserved. In the eigenstate problem, we 
have Lz = Tim, where m is an integer, called the ro- 
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FIG. 7: Effect of the power P2 of the blue-detuned light on 
the total potential CAot . The power of the red-detuned light is 
fixed at Pi = 30 mW. All the other parameters are the same 
as for Fig. 2] 
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FIG. 8: Bound states for the first six levels (n = 0, 1, 2, 3, 4, 
and 5) of the radial motion of a cesium atom in the effective 
potential UcS = il7tot + ?i^(m^ - 1/4) /2Mr^. The rotational 
quantum number is m = 0. All the other parameters are the 
same as for Fig.|l] 



tational quantum number. The centrifugal potential of 
the atom is given by Ud = Ti^{m^ - l/4)/2Mr'^. The 
radial motion of the atom can be treated as the one- 
dimensional motion of a particle in the effective potential 
C^cff = f^tot + ^^cf ■ In Fig. IHl we plot the eigenfunctions 
for the first six levels of the radial motion of the atom 
in the effective potential UcS with the rotational quan- 
tum number m — Q. The energy of the ground state 
is iJo = —2.872 mK. The spacing between the the en- 
ergies of the ground state and the first excited state is 
roughly 23 /iK. The characteristic size of the ground state 
is Ar = 8.8 nm. 



III. TWO-COLOR TRAP WITH LINEARLY 
POLARIZED INPUT LIGHTS 

The formation of a set of trapping minima as a ring 
in the fiber transverse plane (or as a shell in the three- 
dimensional space), shown in the previous section, is due 
to the cylindrical symmetry of all of the components of 
the total potential. In order to confine atoms to vicinities 
of single local points in the fiber transverse plane, we need 
to break the symmetry. An asymmetric optical potential 
can be obtained when one or both trapping light fields are 
linearly polarized at the input. Indeed, it has been shown 
that, due to the thin thickness of the fiber and the high 
contrast between the refractive indices of the silica core 
and the vacuum clad, the intensity distribution of the 
field in a fundamental mode with quasi-linear polariza- 
tion strongly depends on the azimuthal angle, especially 
in the vicinity of the fiber surface [l8l |. 

Regarding the properties of the trapping potential, 
there is no substantial qualitative difference between the 
case where both input fields are linearly polarized and 
the case where one is linearly polarized and the other 
one is circularly polarized. Therefore, we consider here, 
as an example, the case where both input fields are lin- 
early polarized, along the same direction, namely, the x 
direction. In this case, the polarization of each light field 
propagating along the fiber is quasi-linear and the spa- 
tial distribution of the field intensity is not cylindrically 
symmetric Outside the fiber, the time-averaged in- 
tensity of the electric field in mode i is given by pLSL [20j 

\E,\^ = £f{K^iq,r)+w,Kfiq,r)+f,Ki{q,r) 

+ [wiKfiQir) + £.iKo{q,r)K2{qir)] cos2(p}. 

(22) 

Here = 2(1 + Si)/(1 — Si). We use the above field in- 
tensity distributions and Eq. Q to calculate the optical 
potentials Ui and U2 of an atom in the red- and blue- 
detuned fields. The net optical potential J7 = ?7i + [/2 is 
then found to be 

C/(r,(^) = G2{Kl[q2r) + W2Kl{q2r) + hKl{q2r) 

+ [^2^^1(92?-) +6-f^o('?27-)i^2(g2?')]cOs2(^} 

- GiiKliq^r) + w^Kliq^r) + hKl{q^r) 
+ [wiKl{qir) + £,iKo{qir)K2{qir)]cos2ip}. 

(23) 

Unlike the potential (|3J| , the potential H23|) contains addi- 
tional terms that vary with the azimuthal angle Lp. Due 
to these terms, the potential H23|) is not cylindrically 
symmetric. Meanwhile, when the intensity of the blue- 
detuned field is large enough, the potential (|23l) is, like 
the potential ll^Jl, repulsive at short range and attractive 
at long range. For a fixed angle the potential U{r, (p) as 
a function of r must have a local minimum value Umi^p) 
achieved at a point rmi^p) where the repulsive and at- 
tractive forces cancel each other. Due the cylindrical 
asymmetry of the potential, Umif) is not constant in if. 
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FIG. 9: Transverse-plane profile of the total potential C/tot 
produced by the lights that are linearly polarized at the input. 
The hght powers are Pi = 30 mW and P2 = 35 mW. All the 
other parameters are the same as for Fig. |1| 



B 





(a) 


(b) 


4- 




y = 


x = 


2- 




- 




0- 










-2- 










-4- 










C 


) 


1 2 C 


) 


1 2 



X [|J,m] 



Consequently, Um{^) must have a global minimum value 



min Um (if) , which is achieved at a finite number 



of angles ipm. Thus the potential U{r,(p) must have a 
global minimum value Ug achieved at a finite number of 
single points {rmi'fim), 'fim) m the fiber transverse plane. 
Due to the axial symmetry, we have U{r, ip) = U{r, ir+ip). 
Therefore, there must be two minimum points for U{r, (p) 
in the transverse plane. Since the atom can move freely 
along the direction z of the fiber axis, the two minimum 
points in the transverse plane in fact correspond to two 
straight lines for the minima of the trapping potential in 
the three-dimensional space. These potential minimum 
lines can be used to trap and guide atoms along the fiber. 

We use Eq. as well as Eqs. (UHl and (EU to cal- 
culate the total potential C/tot of a ground-state cesium 
atom outside a thin cylindrical silica fiber with two lin- 
early polarized input light fields. We plot in Fig. O the 
spatial profile of J7tot in the fiber transverse plane. The 
figure shows that Utot varies substantially with the az- 
imuthal angle ip, especially in the vicinity of the fiber 
surface. The minimum points in the radial dependence 
of Utot form a ring surrounding the fiber. Unlike the case 
of circular polarization, see Fig. [S] the ring at the bot- 
tom of Fig. 1^1 is not flat with respect to (p. This ring has 
two deepest points, located symmetrically on the x axis 
(at ip = and cp = tt). These points correspond to two 
straight lines, parallel to the fiber axis, for the trapping 
minima in the three-dimensional space. 

To get a closer look at the spatial distributions of the 
potential along different radial directions, we replot in 
Fig. ^1 the potential Utot of Fig. |21 as a function of x at 
y = (i.e., as a function oi r a,t ip = 0) and as a function 
of y at X = (i.e., as a fimction of r at (/s = 7r/2). 
Figures Unta) andHHIb) show that, for Pi = 30 mW and 
P2 = 35 mW, the depths of the potential minima for 
the X and y directions are about —3.2 mK and —2 mK, 
respectively. Thus the depth of the potential minimum 
for the X direction is deeper than that for the y direction. 
Consequently, the local minimum for the x direction is 
the local minimum for the whole transverse plane. The 
difference between the depths of the potential minima 
for the X and y directions is about 1.2 mK. Note that the 



FIG. 10: Total potential fAot as a function of x at y = (a) 
and as a function of y at x = (b). All the parameters are 
the same as for Fig. |2l 
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FIG. 11: Azimuthal dependence of the total potential fAot 
for r = 0.4 um. All the parameters are the same as for Fig.|5] 



radial distances Xm = fmiO) and y„i = r„j(7r/2) of the 
local minima for the x and y directions are almost the 
same, namely, rm{Lp) — 0.4 /xm. 

To get insight into the variations of Utot along the bot- 
tom ring r = rm(sp) = 0.4 /im, we plot this potential 
in Fig. ^] as a function of ipj for a fixed radial distance 
r — 0.4 /im. As seen, along the bottom ring, Utot achieves 
its minimum value at two angles, ip = Q and = tt. The 
difference between the maximum and minimum values of 
Utot along the bottom ring is about 1.2 mK. 

We estimate some important parameters of the trap 
at the minimum point (r = r,„, ip = 0) or [r — r„i,ip — 
tt) in the situation of Fig. |51 We find that the rates of 
scattering due to the light fields at the trapping minimum 

are F^"''^ ^ 25.87 s'^ and F^"''^ = 10.49 s'^. Accordingly, 
the net scattering rate and the coherence time are Tsc — 
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36.36 and Tcoh — 27.5 ms, respectively. The trapping 
lifetime due to recoil heating is estimated to be Ttiap = 
500 s. When we fit the bottom of the potential around the 
minimum point by a 2D harmonic potential, we find that 
the transverse oscillation frequencies are cj^/Stt = 533 
KHz and uJipjl-K = 156 KHz. These frequencies provide 
a ground-state localization of Ir = 8.4 nm and ^ 15.6 
nm. 



IV. CONCLUSIONS 

We have demonstrated that a subwavelength- diameter 
optical fiber carrying a red-detuned light and a blue- 
detuned light can be used to trap and guide atoms out- 
side the fiber. We have shown that, when both the input 
light fields are of circular polarization, a set of trapping 
minima of the potential in the transverse plane is formed 
as a ring around the fiber. In this case, the atoms can 
be confined to a cylindrical shell around the fiber. When 



one or both of the input light fields are of linear polar- 
ization, the potential achieves its minimum value at two 
single points in the transverse plane. In this case, the 
atoms can be confined along two straight lines parallel 
to the fiber axis. Due to the thin thickness of the fiber, 
we can use far-off-resonance fields with substantially dif- 
fering evanescent decay lengths to produce a net poten- 
tial with a large depth, a large coherence time, and a 
large trap lifetime. For example, a 0.2-/xm-radius silica 
fiber carrying 30 mW of 1. 06- /xm- wavelength light and 
29 mW of 700- nm- wavelength light, both fields are cir- 
cularly polarized at the input, gives for cesium atoms a 
trap depth of 2.9 mK, a coherence time of 32 ms, and a 
recoil-heating-limited trap lifetime of 541 s. 
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